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[57] ABSTRACT 

The resolution of an imaging system is greatly enhanced 
by radiating an object with a plane wave field from a 
coherent source variable in either frequency, angle or 
distance from the object, detecting the wave field trans- 
mitted through, or reflected from, the object at some 
point on the image of the object, with or without het- 
erodyne detection, and with or without a lens system. 
The heterodyne detected output of the detector is pro- 
cessed to obtain the Fourier transform as a function of 
the variable for a direct measurement of the amplitude 
and surface height structure of the object within a reso- 
lution cell centered at the corresponding point on the 
object. In the case of no heterodyne detection, only 
intensity data is obtained for a Fourier spectrum. 

30 Claims, 2 Drawing Figures 
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example, a source of ultrasonic waves, and that the 
SUPER-RESOLUTION IMAGING SYSTEM frequency may be held constant while varying the an- 

gle. With off axis illumination in the X direction, the 
ORIGIN OF INVENTION Fourier transform of the wave field with respect to the 


The invention described herein was made in the per- 5 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public law 85-568 
(72 Stat. 435, 42 USC 2457). 
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BACKGROUND OF THE INVENTION 

This invention relates to imaging systems, and more 
particularly to a method and apparatus for greatly im- 
proving the resolution of an imaging system. 

High resolution for imaging systems is a problem due 15 
to the fact that the high frequency detail of the object 
diffracts the incident illumination to angles much 
greater than what the collecting aperture can gather. 
Because of that a considerable amount of information is 
lost in most imaging systems. Therefore, the finite size 20 
of the primary collector presents a fundamental limita- 
tion on the resolving power of any practical imaging 
system. Since the cost and the practical difficulties in- 
volved with building large aperture imaging systems 
increases with the size of this aperture, it has in the past 25 
been of great interest to find alternative means of ob- 
taining a higher resolution with a small aperture. Since 
the case of an infinitely small collecting aperture is 
equivalent to the case of not having any imaging com- 
ponents, imaging or scanning of an object in such a case 30 
is also of great interest. 

SUMMARY OF THE INVENTION 

In a preferred embodiment of the invention, an object 
is radiated off axis (with or without a lens system) by a 35 
plane wave field from a coherent source at an angle, 6, 
and the wave field, or intensity, of the wave field trans- 
mitted (or reflected) by the object is detected at some 
point in the image plane of the object while the spatial 
phase relationship of the radiating wave field is varied, 40 
either by varying the frequency, the angle or the dis- 
tance of the object from the source. The wave field 
(amplitude and phase) or intensity (amplitude) detected 
is processed as a function of the variable to obtain the 
Fourier transform of the wave field to provide an en- 45 
hanced image of the object itself along a line in the 
plane of the illuminating angle. Assuming that the an- 
gle, 0, is in the X-Z plane, where the Z axis is in the 
direction to the object, the image thus produced is along 
the X axis. The position of the object relative to the 50 
imaging system can be moved in a direction normal to 
the plane of the angle, i.e, in the Y axis direction, and 
the process repeated for each new position to obtain a 
two dimensional image. In some applications it is easier, 
or perhaps only possible, to detect just the intensity of 55 
the wave field. Since the intensity is simply the absolute 
square of the wave field, the Fourier transform of the 
intensity with respect to the variable is the autocorrela- 
tion of the image wave field. Thus, while the Fourier 
transform of the intensity does not directly provide the 60 
super-resolution amplitude and phase structure of the 
imaged object, the autocorrelation function provides 
sufficient information to warrant direct measurement of 
intensity. 

A description of the preferred embodiments will em- 65 
phasize the case of a coherent light source, a fixed illum- 
ination angle, 6 , and a variable frequency, but it should 
be understood that the coherent source may be, for 


frequency (inverse wavelength) at a far field region 
point of an image plane yields an image of the object in 
the X dimension with a resolution factor related to the 
span of the frequency variation. Similarly, with a fixed 
wavelength and a variable angle of illumination, the 
Fourier transform of the wave field detected with re- 
spect to the angle of incidence yields an image of the 
object in the dimension of the angle variation with a 
resolution factor related to the span of the angle varia- 
tion possible. When it is possible to use a combination of 
both angular and wavelength variation, one angle in the 
plane of one axis is set to be very large (e.g., $ x =30°) 
while another angle in the plane of another orthogonal 
axis is varied over a small range (e.g., 9 y ±A9 y ). At the 
same time the frequency of coherent source is varied. A 
two-dimensionsl Fourier transform with respect to the 
frequency and angle variations yields an image of the 
object in two dimensions with a resolution factor re- 
lated to the bandwidth of angular and frequency varia- 
tion. It is also possible to use a combination of the angle 
variation in the two directions 6 X , 6 y to provide a two- 
dimensional image. Another possibility is to vary the 
distance to the object along the Z axis instead of vary- 
ing the frequency. 

As used herein, the term “wave field” refers to both 
the amplitude and phase of the waves radiating from the 
source. In the case of electromagnetic radiation, such as 
for example radar. X-ray, or light as in the preferred 
embodiment to be described, the term “electric field” is 
equivalent to the generic term “wave field,” and refers 
to both the phase and intensity of the radiation. The 
generic term wave field is intended to apply to other 
forms of coherent wave propagation, such as ultrasonic 
waves, for example. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of a system implemented to 
carry out the present invention in different modes of 
operation. 

FIG. 2 is a diagram useful in understanding the prin- 
ciples of the invention. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Referring now to the drawings, FIG. 1 illustrates the 
organization of a general, or universal, system for imag- 
ing an object 10 (with an optional lens Li) onto a plane 
11. The system is universal in that it provides a coherent 
source of radiation 12, such as a laser, of tunable wave- 
length under the control of a unit 13, and deflectors 14 
and 15, such as rotatable prisms or electro-optic effect 
devices, for deflecting the illuminating beam in the 
respective X-Z and Y-Z planes under control of a unit 
16. The Z axis is the axis of the beam in the direction of 
the object and the X axis is assumed to be in the plane of 
the drawings. The Y axis is then perpendicular to the 
plane of the paper. Imaging with super-resolution is 
thus possbile with either a variable radiation wave- 
length (frequency) or a variable angle of incidence of 
the radiation transmitted through the object 10, or with 
both a variable wavelength and angle, in which case 
one of the deflectors is set at a large angle (e.g., 30°) 
while the other is varied over a small angle of about 1°. 
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An optional Z-position translator 17 is also provided 
to effectively phase modulate light transmitted through 
the object by varying the position of the object along 
the Z axis instead of varying the wavelength of the 
source 12. Modulation by varying the position of the 
object in the Z plane is fully equivalent to varying the 
angle or the wavelength at normal incidence, as may be 
better appreciated by reference to the diagram of FIG. 
2 which shows the incident coherent light transmitted 


tain the intensity of the radiation source that the differ- 
ence between the feedback signal and a constant refer- 
ence signal (voltage) is maintained at or very near zero. 
However, it is preferable to omit the detector 26 and 
differential amplifier 27, and to merely sample the radia- 
tion just in front of the object using a beam splitter 28 
and a detector 29. The intensity of the radiation thus 
measured is then used in the Fourier transform proces- 
sor 21 to normalize the measurements made through the 
through the object 10 at an angle 0 X to the image plane 10 detector 20. Any variation is intensity caused by the 
11 without the optional lens. FIG. 2 may also be used to deflectors 14 and 15 will thus be compensated. How- 

better understand that when the direction of the inci- ever, it may be advantageous to provide both feedback 

dent light is normal to the object, the transmitted light control of the source and compensation in the proces- 

must then be detected at the image plane off axis for SQr 

super resolution. . lin 15 With or without feedback control of the source, it is 

The electric field at an image point P for the object 10 desirabie to employ the beam splitter 25 to the 

radiation and transmit the sampled radiation to the de- 
tector 20 for heterodyne detection; alternatively an 
independent local oscillator synchronized with the light 
source for heterodyne detection. A beam splitter 30 
may be employed to combine the sampled radiation 


is the electric field transmitted (or reflected) by the 
object over the area of the corresponding resolution cell 
(determined by the resolving power of the imaging 
system) multiplied by the impulse response of the imag- 20 
ing system and integrated over this area. However, the 
electric field exiting the object depends upon the wave- 


length and angle of incidence as well as the phase ™ th ^tion from the object. Reflectors 31 and 32 
(height in the Z axis) and amplitude (length in the X direct the sampled radiation to the beam splitter 30. In 

axis) structure of the object. Thus the value of the elec- 25 the case of feedback control, the reflector 31 would be 
trie field at the image plane 11 varies with the wave- a beam splitter. 

length or angle or incidence on the object. Before describing the operation of the system of FIG. 

Since moving the object (or the source) along the Z I in different modes, the Fourier transform carried out 

axis will vary the phase relationships of the incident in the processor 21 for high-resolution imaging will be 

light as detected at an off axis, far-field point, as can be 30 described. Consider the geometry of FIG. 2 without a 
seen from the diagram of FIG. 2, the value of the elec- lens. A plane wave of a given wavelength is incident on 

trie field at the image point varies with the position of the subject 10 at angles 0 X , 0y with respect to the direc- 

the object along the Z axis. Consequently, the Fourier tion to the object (Z axis). The amplitude transmittance 

transform of this variation provides directly the ampli- of the object is assumed to be a(x,y) and the transmitted 

tude and phase structure of the object within the resolu- 35 radiation is detected at some point P in the far field 
tion cell corresponding to the image point of interest. region or the image plane of the object. Assuming a unit 

It will, in most cases, be more convenient to vary the amplitude illumination, and using scalar diffraction the- 

wavelength (represented by the space between wave- ory in the following calculations, consider only the case 

fronts in FIG. 2) so that it is only necessary to set the of linear polarization. 

deflectors 14 and 15, one for the desired offset angle of 40 The electric field, E n at the output plane of the object 

incidence (e.g., 8 X ), and the other to zero, or to set the io is given by, 
wavelength of the source and vary the angle of inci- 
dence. Consequently, the Z position translator 17 can be p la = a (x^ v )Akcos$xcos 0 yh(xy}^k(xsin 9 x + ystnOy) m 

omitted, as well as one of the angle deflectors 14 and 15, 

except for two dimensional imaging, in which case both 45 where h(x y) is the surface phase transmktance profl i e 


deflectors are required. However, it is desirable to be 
able to move the object in order to image different 
areas. An X-Y position translator 18 is provided for that 
purpose. 

It is also convenient to detect the electric field at 50 
different points in the image plane 11. An X-Y position 
translator 19 is provided to move a detector 2Q in the 
image plane. All position translators provide a Fourier 
transform processor 21 with signals which uniquely 
identify the positions of the object 10 and detector 20, as 55 
well as the values of the variables, namely wavelength, 
and angles 0 X and 8 y in order for the Fourier transform 
to be computed by the processor 21, squared by a 
squarer 22, and the result displayed as an image in a unit 
23 and/or recorded in a unit 24. Position data is re- 60 
corded to provide stereo viewing of the object. 

Since any variations in the intensity of the radiation 
source will result in variations in the electric field mea- 
surements made through the detector 20, it is necessary 
to maintain the intensity of the radiation substantially 65 
constant. That can be done by sampling the radiation 
using a beam splitter 25 and a detector 26 to provide a 
feedback signal to a differential amplifier 27 to so main- 


of the object, and where the time dependent e iat term 
has been suppressed for monochromatic illumination. 

In the super-resolution imaging of an object, three 
possible arrangements are of interest. When an imaging 
system of pupil aperture function P(x,y) is used and the 
impulse response is ^(x,y), the image electric field is 
(for unit magnification) 




(2a) 


where ^(x,y) = JT J\x,y)e. z ^ tox + m ^dx^yo, and 

— QO 

z=the distance of the limiting aperture P(x,y) from the 
image, and * the convolution operation. 

For the case when the Fresnel zone electric field 
scattered by the object is detected at distance z from the 
object, we have 


£f«(Wo) : 


l 

thz 


e ikz e (ik/2z)(X{2 + yit) 


(2b) 

ff + yt) 



5 

-continued 

e (-ik/zKxxo + yyo)j xc jy 
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wavenumber ko- Substituting (4) into (5), obtain the 
result 


eA0,0,7x) = Jff -2^L e /Hxsin»x - yx + z) rect dxdy 

First evaluate Equation (6) with respect to k. Since the 
When the Fraunhofer field is detected at distance z , Q limits of the intergration have been set over the region 
from the object, k<,±Ak/2 by the rect [(k — k<,)/Ak] in the integral, and 


Ef{xo,y 0 ) s 


c ikz e (ik/2z)(xo2 + yo2) 


EJx,y)e ~r~ (xxo + wo) dxdy 


When the field is detected at a single point, the Fraun- 
hofer case is equivalent to the case when the image 
pupil function P(x,y) is a delta function and the imaging 
system has an infinitely poor resolution. The Fresnel 2 q 
case is equivalent to the Fraunhofer case except for the 
fact that the point- at which the field is detected is so 
close to the object that the quadratic phase factors 
across the object become important. In the sub-resolu- 
tion imaging, these quadratic phase factors only con- 2 s 
tribute to a stretching of the image and their effect on 
the actual imaging process can be ignored. Thus only 
the sub-resolution imaging case for the Fraunhofer field 
will be considered and the equivalent calculations for 
the imaging case or the Fresnel field case can be made 30 
correspondingly. Since the system is spatially symmet- 
ric, interchanging (6 x x) with (0 y ,y) does not change the 
essential results. Thus in calculations where 6 y = o, the 
corresponding effect for 0 x =o, 0 y ±o can be easily 
determined. Also for most surfaces, the effect of h(x,y) 35 
is considered equal to zero and will only contribute a 
small random background noise which may be ignored. 

Substituting Equation (1) into Equation (2c), the re- 
sult is 


since the phase term 

e ik(xsin8 x + yx - z) 

is varying at a much greater rate than the 1 A term, the 
variation in the 1A term is ignored in evaluating the 
integral over k. This gives the result 



where sine (x) is defined as sin irx/irx. 

Note that Equation (7) is in the form of a convolution 
and can be written explicitly as 



«W,) = e "S’ ^ r ° 2 + yc2) J f o(x,y) e ikxs ‘ n6x - T- + Wo) dxdy 

— 00 

For the particular point P, at (0,0) on the plane II, Equa- 45 Equation (8) expresses the Fourier transform of the 
tion (3) reduces to electric field with respect to the wavenumber in the far 

field region of an object and shows that it is the electric 
m ( 4 ^ field transmitted by the object convolved with a resolu- 

Ef(Q,0) = jf a(x,y) e ikxsinix dxdy. tion factor related to the bandwidth of the frequency 

1 -co 50 variation available. The width of this resolution factor 

in the x-dimension is 2ir/Ak sin 0 X where 0 X is the angle 
Now evaluate the Fourier transform of Ef(0,0) with of illumination employed. If the wavelength span of 

respect to the wavenumber k. In general, the tunability 1000 A is employed, the illumination angle taken to be 

for the wavelength of the illumination is over some 60°, and the center wavelength taken to be 5000 A, we 

certain region of Ak wavenumbers band-width and at 55 obtain the resolution cell size in the x dimension as 2.9 
some center frequency f 0 which has the corresponding p,m. In the y-dimension the resolution is very poor, 

wavenumber k<,. Thus, in evaluating the Fourier trans- since the image electric field is the integrated field exit- 

form of Ef(0,0) as in (4), take into account this limited ing the object plane over the y-axis. 

bandwidth available. This Fourier transform may, It is possible to proceed to perform a similar calcula- 
therefore, be written as 60 tion for the Fourier transform of the field detected at P 

with respect to the angle of illumination when the 
(5) wavelength of illumination is not varied but the incident 
eKO,0,yx) = f Ef( 0 , 0 )rect ( * 7 .*° ] , angle e x is changed by the span A0* around 0 xo . The 

” oo \ Ak / 2,r pupil function describing this angular bandwidth em- 

65 ployed is assumed to be rect [(sin 6 X — sin 0 XO )/ A(sin 0*)] 
where y^is the image coordinate in the x dimension and where A(sin 0 x )/2=z sin (0*+ A0 x /2)— sin 6 X . The Fou- 

E/fO.O) has been producted by rect [(k-K<,)/Ak] to take rier transform of Ef(0,0) is defined, using the same 

into account the limited bandwidth available around the notation used in Equation (6) as follows 
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Now taking the two-dimensional transform of this 
ef ( 0 fi,p x ) = ( 9 ) function as defined by 

/ sin0* — sin9„ \ , . . - . 

/ EM0)rect j e~**»*&* d(sin8 x ) eF (y x ,p y ) = 


Substituting Equation (4) into Equation (9), evaluating, •'T. «, ' ' V 

and simplifying the result, the function ef(o,o,p x ) is , v 

found to be given by the convolution e-®w - a*W*ji>py d f — j (/(sine,,) 

in \ / 


( 13 ) 

/•/.“ ( k — k 0 \ ( sinSj, \ 

// M j rec \ ^( S ineJ J 


e/fo.o. Y' ) = 


a(*y)» 


giifcz — /2?r sindxo 
ikz 


^ I0 ' 1 and substituting Equation (12) into Equation (13) the 

[ _ 1 following equation is obtained: 

A(sin0Jr) —■ I (Asm#,) 

ep(y x ,py) = fill -gfcjfl. . e /*(xswi0* -f — yx + z) — il-n{sindy)py 

— CO 


(14) 


rec ' ( lL nr ) rect ( “aSt ) rf ("5r) d ^ nd y> dxd y 


Evaluating Equation (14) over only the k variable, ne- 
glecting the 1/X variation compared to the phase term 
over the region k 0 ±Ak, yields 


efiyx.Py) = Iff q (^Y? e iko(xsinO x + ysinOy — yx + z) — i27r(sindy)p y 


( 15 ) 


( Ait 

^ V sinO, ^ . F 

^ 2t r 

J rcC 'l, A(sin0,) )* mc [ 


The Fourier transform of the electric field detected at 
point P with respect to the angle of incidence is thus 
obtained as the convolution of the electric field exiting 
the object plane I with the function {(eikz/i\z^ e —aw(x/\) 
sin 0 X °( A Sin 6x ) sine [A(sin 0^)x/X]} which has the width 
given by X/A(sin B x ). Thus, the resolution of the image 
obtained in this fashion, for an angular width of ±5° 
centered around normal incidence, and for an incidence 
wavelength of 5000 A, is 2.9 pm in the x-dimension. 
Again, the resolution in the y-dimension is very poor 
since the field is integrated over the whole y-axis. 

It is possible to obtain a high resolution in both di- 
mensions by using a combination of both wavelength 
and angular variation monitoring. Consider the case 
when the illumination angle in the direction of the x axis 
is 0 X , and the angle that the y component of the k vector 
for the illumination makes with the normal to the object 
is By, and assume that B x > >B y . In computing the Fou- 
rier transform of the electric field with respect to the 
variation in the wavelength of the illumination and in 
the angle 0 y , the electric field exiting plane I is given by 


EJ,x,y)=a(x,y)e> kx « x+ikysin By. 


(U) 


35 


and the corresponding electric field detected at point P 
is written as 


ikz 00 ( 12 ) 

E *w- ir If a(x,y)e ikxsin ^ x + ‘kysinty j x( / y 


Recording the value of E/70,0) while the incident 
wavevector is varied from ko — Ak/2 to ko+Ak/2 and 
the angle of incidence in the y direction from — AB y /2 
to +AB y /2. yields the function E/^0,0) rect 
[(k— k»)/Ak] rect [sin 0 > /A(sin 0 y )]. 


60 


65 


Since it is assumed that Byo=0, and that the variation of 
By is so small that A B y < < B x , it is possible to ignore y sin 
By compared to the x sin 6 X term. The sine 
[Ak/27r(yjc— x sin B x — y sin B y —z)] can be approxi- 
mated by sine [Ak/27r(y x — x sin 6 x —z)]. Using this 
approximation to evaluate Equation (15) yields 


40 


eKy x ,py) = 


( 16 ) 


45 


JJ dxdy A(sin0j,) J eW**" - yx + z) 

sinc ["fjr" ~ xsin6x ~ z ) J sinc ^X(sinflj,) ^p y - 


which can be rewritten as the convolution 


50 


^sinS* = a(x,yr e -iMxsin0x + z) A(sin ^) J 

sinc ^ (xsin6 x - z) J sinc j^A(sinO^) — J . 


55 


Thus, provided By+A9 y <. <B X , the resulting two-di- 
mensional Fourier transform with respect to the wave- 
number k and the angle of illumination variation B y , 
would give the electric field transmitted by the object 
10 convolved with the impulse response factors for the 
x and y dimensions, the widths of which are 1/Ak sin B x 
and X/A(sin By) in the x and y dimensions respectively. 
Thus, as in previous considerations, for X=5000 A, 
AX=1000 A, B x = 60°, B y = 0°±5° this would corre- 
spond to a resolution of 2.9 pm in the x and y dimen- 
sions respectively. Consider now the scanning applica- 
tion where high resolution is obtained in one dimension 
and the object translated with respect to the illumina- 
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tion in the other direction. In order to do this, assume 
that the illumination spot on the object is given by the 


when 6 x =8y= 0, h(x,y)±0, and the integral (5) is evalu- 
ated. The result is 


CO 

ep(o,o,y x ) = Jf o(jc,y)e - '*°b' ~ 2.h(x,y) - z ) sinc[dt/2ir(y^ — 2h(x,y) — z)] 


function rect (x/2L*) rect (y—a/2L y ) which is a rectan- 
gle of width 2L*, 2h y 'm the x and y dimensions respec- 
tively and centered at (0, a). While the wavelength or 
angular variation is used to obtain the high resolution 
along the x dimension the object is translated with re- 
spect to the illumination in the y direction, i.e., a is 
varied. For any particular value of a, use (7) to provide 
the Fourier transform with respect to the wavenumber 
of the electric field detected at the point P. This is given 
by, for the object transmittance now given by 


For Ak large, ef(y x ) has non-zero values only for 
y=2h. This point may be clarified further by assuming 
10 that there is only one phase scatterer on the object and 
the function h(x) is described by h(x)=h 0 rect(x— x- 
o)/w 0 where (x 0 ,o) is the position of this scatterer of 
height h 0 and width w„. 

The function efiy x ) will have a non-negligible value 
15 only at y=2h 0 and the magnitude of this value will be 


( a(x,y)e i,cxsin0x rect ^rect 

co a + Ly / \ MxtinO* - y + z)( a l \ 

eftO.O.yj) = //^ a _ l > [— J (tSt J 

sine (y x ~ 2 - xsinff x ) ~^dxdy. 


Differentiating (18) with respect to a and rearranging 
terms, we have the result, 30 


a{x, a)rect 


( )*' 


-iko(sin9x + z) 
fhz 


•(^-)sinc - 


• (z + xsinfl*) 


de F ( x \ 

-fa- (p,o.y x ) + a(x,a - 2L y )rect I -jj— J 

(■^) sinc ( z + xsin Qx) Js(a) 


* e iko(sin$x + z) 


From Equation (19) it is evident that, provided a(x,y) is 
known for the region 0<y<2L^, the function a(x,y) 
may be determined from Equation (21) for all values of 
y for any given L y . Thus, for scanning applications, ^5 
while high resolution may be obtained in one dimension 
by varying the wavelength or angle of illumination, the 
other dimension may be resolved if the object is trans- 
lated with respect to this axis and if the object transmit- 
tance is known in the initial spot size region. 30 

The electric field at the point 03,0) also depends upon 
the distance z of the detector plane from the object. 
Evaluating the integral 


* l/z - 1/Z„ 

j Ep(J$,o) reel &(i/z) " 


-ajri)(i/z) d(i/z), 


for z 0 a mean distance, E/03,0) as defined in Equation 
(3) and 8 X = 0, /3/z> > 1 60 


/ a{x)rect( -■ — *° - )dx. 
— co 0 


For a multitude of scatterers of height h r , width w r on 
the object, the function e/f'/x) will be a series of peaks 
at yx=2A r and the heights of these peaks will be equal to 


/ a(x)rect( )dx. 

For a continuous function describing h(x) the function 
SFiy x ) will be some continuous function for various 
values of y x depending on the fraction of the total func- 
tion that has the height equal to y x produced by the 
values of a(x,y) at those particular y=2h. By placing a 
mask, with step function heights of equal widths and 
different values placed at different spatial points, in 
contact with the object a(x,y), the values of a(x,y) can 
be sampled by measuring the Fourier transform of the 


rnfik -ik/zdxP - 

— ) = <x,yy 


■ (A( 1 /z))sinc[A( 1 /z)/3/\x - P/2)] 


While Equations (3) to (20) have been evaluated with 
the assumption that the spatial effect of the object phase 
transmittance h(x,y) is zero, consider now the case 


electric field detected versus wavelength and noting 
that the position of the peaks on this transform corre- 
spond to the heights on the phase mask and the value of 
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these peaks to the value of a(x,y) at the corresponding 
points on the mask. 

While Equations (8), (10), (17) and (19) have been 
calculated for the case when point P is at (0,0), it is 
interesting to note the effect of making the imaging 
measurement at some point other than the origin. To 
illustrate, consider- the derivation of Equations (7) and 
(10) using the values x 0 = — /? and y o =0 in (3). Proceed- 
ing with similar steps leading from Equations (3) to (8), 
it is found that the Fourier transform with respect to the 
wavenumber of the field detected at (— /J,0) is given by 
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nately the object a(x,y) can be written as a matrix in 
which case m is a matrix and m _1 is an inverse matrix. 

In the practical implementation of the theory in- 
volved in the analysis here, a coherent illumination 
5 tunable in wavelength or angular direction and a 
method for detecting the far field electrical field scat- 
tered by the object, such as by a heterodyne measure- 
ment, would be sufficient to provide the necessary in- 
formation to produce a high resolution image of the 
10 object of interest. For optical wavelengths, an electro- 
optically tunable dye laser or a laser in combination 
with an acoustic-optic or electro-optic deflector may be 


e/^0,0, ^sin0* + -f-^ = a(x,y) • e - '* 0 + "t") + (* + 

(-££-) sinc [-f£- x ( S me, + + ( z + ~fr)] 


e/^0,0, = 


Similarly the equivalent of Equation (10) becomes 

( 23 ) 25 

a(x,y) • e~ ik + lr 

A(sin0 x ) sine ^AfsinO*) j-j- J. ^0 


Note from Equation (22) that detecting the wave- 
length dependent electric field at a point other than the 
origin simply shifts the image by /3 2 /2z with respect to 35 
the image formed from the field at the origin. Also the 
resolution cell size of the image has now changed from 
2rr/k sin 0 X to 27rA(sin0 x -|-/3/z). Thus, even if the il- 
lumination angle is zero, i.e., if sin 0 X =O, detection of 
the electric field would still produce some image, the 40 
resolution of which is proportional to the angle sub- 
tended by the detector, fS/z, to the object. This is re- 
lated to the resolution for a lens of aperture /?. From 
Equation (23), the off center measurement of the elec- 
tric field, while the angle of illumination is changed, 45 
does not perturb the measurement except to change the 
phase factor involved with the impulse response func- 
tion. 

Equation (22) is also applicable, if 0 X =O, when any of 
the standard interferometric techniques to make inco- 50 
herent holograms are used; and field at the point (/3,0) is 
detected with respect to wavelength. In the case when 
a band of frequencies illuminates an object at the same 
time, the electric field at a point on the far-field region 
detected as a function of time is the image electric field 55 
given by Equation (8) since the time and frequency 
variables are Fourier transform variables with respect 
to each other. 

Given for some function a(x,y), the operation 
g(xi,yi)= f f a(x,y)m(x,y;x \,y \ )dxdy, the function a(x- 60 
0 ,y 0 ) may be retrieved from g(xj,yi) by using the opera- 
tion //g(xi,yi)m- 1 (x,y,y JC ,y J ,)dxidyi where the func- 
tion m- 1 (x,y;y JC ,y J ,) is defined by / $m(x,y;x\,yi)m~\x- 
by\;yx,yy)dx\dy}=d(x~y x )d(y-yy). The function m in 
the situations described in the subresolution system 65 
corresponds to e' 2 ’ rxxl and m~ 1 to e-i^yxxi where x; is 
1A when wavelength variation is used for the imaging, 
sin 0 X or sin 0 y when angular variation is used. Alter- 


( 22 ) 


used to produce rapid tuning of the wavelength and 
angle of illumination, while the electric field at point P 
may be detected by using a local oscillator to provide 
the coherent detection scheme. 

The Fourier transform operation to reduce this signal 
to an image may be performed by a digital fast Fourier 
Transform Processor. An alternate technique would be 
to amplitude and phase modulate the output of a laser 
by an electro-optic cell with the signal and detect this, 
using heterodyne detection, with an integrating time 
longer than the signal duration. (The local oscillator 
may be a tunable laser.) The detector output obtained as 
a function of the frequency difference of the local oscil- 
lator and the signal laser would be the desired Fourier 
transform. 

Equations (8), (10) and (17) are the results for the 
Fourier transform of the electric field, with respect to 
the angle or wavelength of illumination incident on 
some object at a point in the far field region of the 
object. They represent the image electric field of the 
object itself, the resolution of this image depending 
upon the angular and wavelength bandwidth through 
which the illumination can be tuned. 

The foregoing discussion of the system of FIG. 1 in 
different modes has been generally in reference to the 
system without the optional lens Li, particularly in 
respect to the geometry of FIG. 2 and the discussion 
relation thereto which followed. For the corresponding 
geometry in the case of imaging with a lens, and a dis- 
cussion of the present invention as it applied to that 
case, see a paper by the present inventor titled Sub- 
Resolution Imaging: Use of Wavelength and Angular 
Diversity Techniques published in Applied Physics 
Letters, Vol. 28, No. 8, at pages 451 to 453 (Apr. 15, 
1976). That paper also illustrates the applicability of the 
invention to applications employing reflected wave 
field detection. 

In summary, the Fourier transform of the wave field 
is obtained by a detector at a point on the image plane 
of an object as the angle, wavelength or distance from 
the object of the coherent wave source is varied to 
provide through an inverse Fourier transform calcu- 
lated by a processor enhanced image data of the object 
which the resolution cell centered at the corresponding 
point on the object. Alternatively, only the intensity of 
the wave field may be detected as a Fourier spectrum at 
a point on the image plane, and the Fourier transform of 
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the intensity is the autocorrelation of the transform of 
the wave field variation with wavelength. Thus, it is 
possible to use the output of the detector without het- 
erodyne detection in order to obtain only intensity 
(wave amplitude) data for a Fourier spectrum of the 5 
image (i.e., the square of intensity of the Fourier trans- 
form that would otherwise be obtained by heterodyne 
detection), or to use the output of the detector with 
heterodyne detection to obtain both the amplitude and 
phase data of the image for a Fourier transform. This is 10 
so because the output of the heterodyne detector is the 
Fourier transform of the object at a point on the image 
plane, so that if the Fourier transform processor is not 
used, the result is a spatial Fourier transform of the 
object; and if the processor is used to obtain the inverse 15 
Fourier transform, a wave field image at a point on the 
image plane of an object is produced. Upon squaring 
that wave field image, an intensity image is produced 
for display, as on a cathode ray tube display unit. In 
other words, the squarer provides the square of inten- 20 
sity in the image wave field, and eliminates the phase 
data. Also a system without a processor may be of inter- 
est in some applications since the detector may provide 
the Fourier transform of the object if both modulation 
and heterodyne detection is used, as described hereinbe- 25 
fore, and the Fourier spectrum if only intensity detec- 
tion is used, i.e., if elements 25, 31, 32 and 30 are omit- 
ted, as just noted. While the preferred embodiments 
described herein have been for the case of detecting 
radiation through the object, it is obvious that they can 30 
be modified for the case of reflection from the object. 
Consequently it is intended that the claims cover such 
modifications and variations. 

What is claimed is: 

1. A method of enhancing the resolution of a system 35 
for imaging an object at a plane comprised of the steps 
of 

radiating said object with a plane wave field at an 
angle with respect to an axis normal to the object 
and in the direction to the object from a coherent 40 
source variable in a parameter selected from fre- 
quency, angle, or distance to the object, to achieve 
phase variable radiation of the object, 
detecting the wave field or intensity of radiation at a 
selected point in said plane as a function of the 45 
variable parameter selected while said selected 
parameter is being varied, and 
computing a Fourier transform of the detected wave 
field or intensity at said selected point with respect 
to the variable parameter. 50 

2. A method as defined in claim 1 wherein said se- 
lected point of the imaging system includes the whole 
object. 

3. A method as defined in claim 1 wherein said se- 
lected point of the imaging system includes a single cell 55 
of the whole image and said cell includes an area of the 
image determined by the resolution of the imaging sys- 
tem. 

4. A method for enhancing the creation of an image at 

a plane comprised of the steps of 60 

radiating an object with a plane wave field from a 
coherent source at a selected angle 0 X with respect 
to a Z axis normal to the object and in the direction 
to the object from said source in a plane defined by 
an X axis, where said X and said Z axes are perpen- 65 
dicular and perpendicular to a Y axis to define an 
X-Y plane parallel to said object and an X-Z plane 
of said angle, 9 X , 
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varying the spatial phase relationship of the radiating 
wave field by varying at least one parameter se- 
lected from the wave field frequency and the se- 
lected angle of radiation, 

detecting the wave field of radiation at said selected 
point as a function of the variable parameter se- 
lected while said selected angle 0 X is being varied, 
and 

computing a Fourier transform of the detected wave 
field or intensity of a selected point in said plane 
with respect to the variable parameter for synthesis 
of a high resolution image at said point. 

5. A method as defined in claim 4 wherein said Fou- 
rier transform is of the image of the whole object along 
a line in the plane of the angle. 

6. A method as defined in claim 4 wherein said Fou- 
rier transform is of a cell of the image, and said cell 
includes an area of the object determined by the resolu- 
tion of the imaging system. 

7. A method as defined in claim 4 wherein both the 
radiation frequency and angle are varied, the further 
limitation of fixing the angle of radiation in the plane of 
one axis normal to an axis in the direction to the object 
and limiting the variation in the angle in the plane of 
another axis normal to said one axis and the axis in the 
direction of the object over a small range while varying 
the frequency, and including the computation of Fou- 
rier transforms with respect to both the frequency and 
angle variations to obtain an image of the object in two 
dimensions. 

8. A method as defined in claim 4 wherein only the 
angle of radiation is varied, first in a plane of one axis 
normal to an axis in the direction of the object, and then 
in a plane of another axis normal to said one axis and the 
axis in the direction of the object to obtain an image of 
the object in two dimensions. 

9. A method as defined in claim 4 wherein the Fourier 
transform computation is achieved by modulating the 
amplitude and phase of radiation from said coherent 
source with a periodic signal and, using heterodyne 
detection with an integrating time longer than the signal 
period, objecting a signal as a function of the frequency 
difference of the coherent source and the local oscilla- 
tor that is the desired Fourier transform. 

10. A method for enhancing the creation of an image 
in an imaging system comprised of the steps of 

radiating an object with a plane wave field from a 
coherent source at a selected angle 0 X with respect 
to a Z axis normal to the object and the direction to 
the object from said source in a plane defined by an 
X axis, where said X and said Z axes are perpendic- 
ular and perpendicular to a Y axis to define an X-Y 
plane parallel to said object and an X-Z plane of 
said angle, 0 X , 

varying the distance of the object from the source in 
order to obtain high frequency components of an 
area at a selected point in the image for Fourier 
synthesis, 

detecting the wave field or intensity at said selected 
point as a function of the variable distance while 
said variable distance is being varied, and 

computing a Fourier transform of the detected wave 
field or intensity of said selected point with respect 
to the variable distance for synthesis of a high reso- 
lution image at said point. 

11. A method as defined in claim 10 wherein said 
Fourier transform is of the whole object along a line in 
the plane of the angle. 
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12. A method as defined in claim 10 wherein said 
Fourier transform is of a cell, and said cell includes an 
area of the object determined by the resolution of the 
imaging system. 

13. Apparatus for enhancing the resolution of a sys- 
tem for creating an object at a plane comprising 

means for radiating an object with a plane wave field 
from a coherent source variable in a. parameter 
selected from frequency and angle of illumination 
to the object to achieve phase variable radiation in 
order to obtain high frequency components of an 
area at a selected point in said plane for Fourier 
synthesis, 

means for detecting the wave field or intensity at said 
selected point as a function of the variable parame- 
ter selected while said selected parameter is varied, 
and 

means for computing a Fourier transform of the de- 
tected electric field or intensity of said selected 
point with respect to the variable parameter. 

14. Apparatus es defined in claim 13 wherein said 
selected point of the imaging system includes the whole 
object. 

15. Apparatus as defined in claim 13 wherein selected 
point of the imaging system includes a single cell of the 
whole image and said cell includes an area of the image 
determined by the resolution of the imaging system. 

16. Apparatus for enhancing the creation of an image 
in an imaging system comprising 

means for radiating an object with radiation from a 
coherent source at a selected angle 0 X with respect 
to a Z axis normal to the object and in the direction 
to the object from said source in a plane defined by 
an X axis, where said X and said Z axes are perpen- 
dicular and perpendicular to a Y axis to define an 
X-Y plane parallel to said object and an X-Z plane 
of said angle, 0 X , 

means for varying the spatial phase relationship of the 
illuminating radiation by varying at least one pa- 
rameter selected from the radiation frequency and 
the selected angle in order to obtain high frequency 
components of an area at a selected point in the 
image for Fourier synthesis, 

means for detecting the wave field or intensity of 
radiation at said selected point as a function of the 
variable parameter selected while said selected 
parameter is varied, and 

means for computing a Fourier transform of the de- 
tected electric field or intensity of said selected 
point with respect to the variable parameter for 
synthesis of a high resolution image at said point. 

17. Apparatus as defined in claim 16 wherein said 
Fourier transform is of the image of the whole object 
along a line in the plane of the angle. 

18. Apparatus as defined in claim 16 wherein said 
Fourier transform is of a cell of the image, and said cell 
includes an area of the image determined by the resolu- 
tion of the imaging system. 

19. Apparatus as defined in claim 16 wherein both the 
radiation frequency and angle are varied, the combina- 
tion including means for fixing the angle in the plane of 
one axis normal to an axis in the direction to the object 
and limiting the variation in the angle in the plane of 
another axis normal to said one axis and the axis in the 
direction of the object over a small range while varying 
the frequency and including means for the computation 
of Fourier transforms with respect to both the fre- 


16 

quency and angle variations to obtain an image of the 
object in two dimensions. 

20. Apparatus as defined in claim 19 wherein only the 
angle is varied, first in a plane of one axis normal to an 

5 axis in the direction of the object, and then in a plane of 
another axis normal to said one axis and the axis in the 
direction of the object to obtain an image of the object 
in two dimensions. 

21. Apparatus for enhancing the creation of an image 
10 in an imaging system comprised of 

means for illuminating an object with radiation from 
a coherent light source at a selected angle 0 X with 
respect to a Z axis normal to the object and in the 
direction to the object from said source in a plane 
15 defined by an X axis, where said X and said Z axes 
are perpendicular and perpendicular to a Y axis to 
define an X-Y plane parallel to said object and an 
X-Z plane of said angle, 0 X , 

means for varying the distance of the object from the 
20 source in order to obtain high frequency compo- 
nents of an area at a selected point for Fourier 
synthesis, 

means for detecting the electric field or intensity of 
illumination at said selected point as a function of 
25 the variable distance while said variable distance is 
varied, and 

means for computing a Fourier transform of the de- 
tected electric field or intensity of said selected 
point with respect to the variable distance. 

30 22. Apparatus as defined in claim 21 wherein said 

Fourier transform is of the image of the whole object 
along a line in the plane of the illuminating angle. 

23. Apparatus as defined in claim 21 wherein said 
Fourier transform is of a cell of the image, and said cell 

35 includes an area of the image determined by the resolu- 
tion of the imaging system. 

24. A method of enhancing the resolution of image 
data comprised of the steps of 

radiating an object with a plane wave field from a 
40 coherent source variable in at least one parameter 
selected from frequency, angle of radiation and 
distance to the object to effectively achieve phase 
variable radiation of the object in order to obtain 
high frequency components of an area at a selected 
45 detection point, and 

detecting the intensity of radiation at said selected 
point as a function of the variable parameter se- 
lected while said selected parameter is varied to 
obtain a spatial Fourier spectrum of the detected 
50 intensity at said selected point with respect to the 
variable parameter for autocorrelation data of the 
image at said point. 

25. A method as defined in claim 24 wherein said 
selected detection point of the imaging system includes 

55 the whole object. 

26. A method as defined in claim 24 wherein said 
selected detection point of the imaging system includes 
a single cell of the whole image and said cell includes an 
area of the image determined by the resolution of the 

60 imaging system. 

27. A method of enhancing the resolution of image 
data comprised of the steps of 

radiating an object with radiation from a coherent 
source variable in at least one parameter selected 
65 from frequency, angle and distance to the object to 
achieve phase variable radiation of the object in 
order to obtain high frequency components of an 
area at a selected point, and 
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detecting the intensity of radiation at said selected 
point as a function of the variable parameter se- 
lected while said selected parameter is varied by 
using a heterodyne detector in order to obtain a 
spatial Fourier transform. 5 

28. A method as defined in claim 27 wherein said 
selected point of the imaging system includes the whole 
object. 

29. A method as defined in claim 27 wherein said 
selected point of the imaging system includes a single 10 
cell of the whole image and said cell includes an area of 
the image determined by the resolution of the imaging 
system. 

30. A method of enhancing the resolution of a point 
on an image plane of an imaging system comprised of IS 
the steps of 
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radiating an object with a wave field variable in a 
. parameter selected to vary the spatial distribution 
of the wave field at a selected region on the object 
corresponding to a point on said image plane, 
detecting the value of the wave field at said corre- 
sponding image point as a function of the variable 
parameter selected while said selected parameter is 
varied, and 

determining the function which, when correlated 
with a wave field detected at said corresponding 
image point as a function of the variable selected, 
results in a response corresponding to a wave im- 
pulse at a position corresponding to the location of 
said point on said object and correlating the wave 
field at the image point with this function. 
***** 
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